powered by
It returns the maximum likelihood estimate of the Matrix Fisher parameter F(3x3).
matrixfisher.mle(X)
The components of \( svd( \bar(X) ) \).
An array containing rotation matrices in SO(3).
Anamul Sajib & Chris Fallaize.
R implementation and documentation: Anamul Sajib <sajibstat@du.ac.bd> & Chris Fallaize.
Prentice M. J. (1986). Orientation statistics without parametric assumptions.Journal of the Royal Statistical Society. Series B: Methodological 48(2).
rmatrixfisher
F <- 10^(-1) * matrix( c(85, 11, 41, 78, 39, 60, 43, 64, 48), ncol = 3 ) ### An arbitrary F matrix X <- rmatrixfisher(5000, F) matrixfisher.mle(X) svd(F)
Run the code above in your browser using DataLab